Frustration and glassiness in spin models with cavity-mediated interactions 
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We show that the effective spin-spin interaction between three-level atoms confined in a multimode 
optical cavity is long-ranged and sign-changing, like the RKKY interaction; therefore, ensembles of 
such atoms subject to frozen-in positional randomness can realize spin systems having disordered 
and frustrated interactions. We argue that, whenever the atoms couple to sufficiently many cavity 
modes, the cavity-mediated interactions give rise to a spin glass. In addition, we show that the 
quantum dynamics of cavity-confined spin systems is that of a Bose-Hubbard model with strongly 
disordered hopping but no on-site disorder; this model exhibits a random-singlet glass phase, ab- 
sent in conventional optical-lattice realizations. We briefly discuss experimental signatures of the 
realizable phases. 



Realizing models of magnetic phenomena and explor- 
ing their phases has been a central objective in ultracold 
atomic physics since the advent of optical lattices [1]. 
Such models (e.g., the Hubbard and Heisenberg mod- 
els) have been of long-standing theoretical interest as 
they are believed to offer minimal descriptions of strongly 
correlated materials [2]. Unlike real materials, ultracold 
atomic systems offer the prospect of realizing the theoret- 
ical models exactly; because many of the models are not 
solvable, it is hoped that ultracold-atomic realizations 
will shed light on their properties. The central effort to 
realize magnetism, to date, has focused on the fermionic 
Hubbard model [3]; however, its magnetic ordering tem- 
perature is too low to be readily achievable in current 
experiments. These difficulties have stimulated an inter- 
est in alternative paths to quantum magnetism [4-6], of 
which the present work is an example. 

We introduce a scheme for realizing magnetism — 
involving A-type three-level atoms [see Fig. 1 (a)] trapped 
in a multimode optical cavity that differs from previous 
schemes in an essential respect, viz. the range and struc- 
ture of interactions. Whereas previous schemes have in- 
volved contact or dipolar interactions [6], the spin-spin 
interactions in our scheme, being mediated by cavity 
modes, are both long-ranged (indeed, infinite-ranged for 
a single-mode cavity) and oscillatory in sign. In these 
respects, they resemble the Ruderman-Kittel-Kasuya- 
Yosida (RKKY) interaction [7], which underlies, e.g., 
the physics of heavy-fermion materials [8] and metal- 
lic spin glasses [9, 10]. The present work is concerned 
chiefly with the latter class of systems, and, in particular, 
with the fact that long-range, sign-changing interactions 
between spins facilitate the realization of various frus- 
trated and bond-disordered models. (Analogous realiza- 
tions have also been proposed, e.g., in photonic band-gap 
systems [11] and Coulomb crystals [12].) 

The elements of our scheme are A-type atoms (i.e., 
atoms with the level structure shown in Fig. 1) dressed 
by a configuration of laser and microwave fields suggested 



in Ref. [13]; the atoms are assumed to be tightly con- 
fined near fixed, random positions inside the cavity. (Al- 
ternatively, the spins might arise because of nitrogen- 
vacancy centers in diamond [14], distributed randomly in- 
side a multimode cavity.) We show that the effective spin 
Hamiltonian is a variant of that studied in Refs. [10, 15- 
17]. By adapting the results of Refs. [10, 16], we show 
that, depending on the number of spins per strongly- 
coupled cavity mode, the low-temperature phase is either 
a spin glass or a superradiant phase [18-20] (analogous to 
a cavity-mediated crystal [21, 22]); we discuss how these 
phases can be distinguished experimentally. In contrast 
with condensed-matter realizations, the systems consid- 
ered here allow one to access both regimes in the same 
system, by changing the mirror spacing and thereby tun- 
ing the number of active cavity modes [23]. We note, 
moreover, that for quantum spins, the effective Hamilto- 
nian can be mapped onto a Bose-Hubbard model possess- 
ing strictly off-diagonal disorder [24]. Unlike the diago- 
nally disordered Bose-Hubbard model [25, 26], the off- 
diagonally disordered version exhibits multiple distinct 
insulating phases, including a Mott glass phase and a 
random-singlet glass phase [24, 27, 28], neither of which 
has been experimentally observed so far. 

Model. We consider A-type atoms whose lower lev- 
els (which will be our two spin states, |+) and |— )) are 
separated by a microwave transition whereas the excited 
level, |e), is separated from both by an optical transition. 
The |±) states are assumed to be tightly confined at the 
intensity extrema of trapping lasers that are far detuned 
from the |±) — > |e) transition; i.e., the atomic positional 
degrees of freedom are assumed to be frozen out. (Thus, 
the physics considered here differs from that of mobile 
spinful atoms [29-31].) Disorder can be introduced using 
diffusers (see, e.g., Ref. [26]). The atoms are confined in 
an optical cavity having multiple degenerate modes, at a 
frequency red-detuned from the |+) — > |e) transition by 
A ~ 1 GHz; other modes are typically farther-detuned 
(e.g., by <~ 15 GHz for a 1 cm cavity). Additionally, the 
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FIG. 1. (a) Level structure of three-level A atoms, dressed by 
a pump laser at frequency ujl, cavity mode(s) at frequency 
u)c, and a microwave field represented by h. The detuning 
from two-photon resonance, 5, is assumed to be much smaller 
than the detuning of laser and cavity photons from the atomic 
transition, A. (b) Proposed experimental setup. Atoms are 
tightly trapped by trapping lasers, which are far detuned from 
the atomic transition, and pumped transversely. Spins are 
self-organized as discussed in the text for a single-mode cav- 
ity, with a sinusoidal mode function as depicted: spins at even 
antinodes interact ferromagnetically with spins at other even 
antinodes, but antiferromagnetically with spins at odd antin- 
odes. Spin-spin interactions are strongest for spins trapped 
at antinodes; therefore, ordering is strongest at antinodes and 
weakest at nodes. 

atoms interact with a pump laser oriented transverse to 
the cavity axis, red-detuned from the |— ) — ¥ \e) tran- 
sition by A + 5, where 5 ~ 10 MHz is the detuning 
from two-photon resonance. (Note that the aforemen- 
tioned setup generalizes to systems possessing manifolds 
of ground and excited atomic states, rather than two 
ground states and one excited state, provided that — as 
in dysprosium [32] — the two manifolds have similar g- 
factors.) The microwave |+) <-> |— ) transition is driven 
at a weak Rabi frequency that, as we shall see, acts as 
an effective magnetic field. 

Under these conditions, the spin-spin interactions can 
be understood as follows: an atom in the |— ) state can 
scatter a laser photon into a cavity mode, thus changing 
its state to | +) ; this virtual cavity photon, being 5 higher 
in energy than laser photons, is reabsorbed into the laser 
after a time ~ 1/5. The reabsorption involves flipping 
the state of a |+) atom (typically a different one from 
the initial atom) to |— ). This entire process generates an 
effective interaction of the form M(xj, Xj)a l + a :l _/(A 2 5) 
between two atoms, where Xj is the position of the i th 
atom and M is a matrix element (derived below) de- 
pending on the cavity mode(s). 

We further assume that the cavity photon leakage rate 
per mode, n <C 5, and also that the atomic-excited-state 
decay rate, 7 -C A. In this "dispersive" regime, the con- 
servative virtual-excitation processes fall off as 1/5 and 
1 /A respectively, whereas the dissipative processes fall off 
as k/5 2 and 7/ A 2 respectively. As argued in Ref. [22], 



the effect of dissipation in this regime can be understood 
in terms of heating, and need not be explicitly included in 
the Hamiltonian. Generally, dissipation does not change 
the mean-field properties even beyond this regime [33]; 
we shall revisit this issue in future work. (Note that a 
weak microwave field must be applied to prevent the spin 
population from being pumped entirely into the |+) state; 
however, this field is comparable in strength to the decay 
processes, and thus much weaker, for small k/5, than the 
interaction terms.) Thus, we neglect, in this work, issues 
such as the nonequilibrium growth of entanglement [13]. 

Hence, upon adiabatic elimination [22, 34] of the state 
|e), the Hamiltonian H of the atom-light system takes 
the form 



H = H ^+*^2 a ^ a ala a + — ^ 5a(xi)c r lal+h.c, (1) 

a, i—l 

where uj a is the frequency of cavity mode destroys a 
cavity photon; f2 is the strength (i.e., Rabi frequency) of 
the pump laser; <7 Q (Xj) describes the coupling to mode a 
at the position Xj of atom i; and the a operators are Pauli 
matrices acting on the atomic ground-state manifold. 
One can rewrite the coupling g a {^i) as <jE a (x,-), where g 
is an overall coupling strength (assumed to be the same 
for all strongly-coupled modes) and 3 a (xj) a normalized 
mode profile. The terms in £f at = J2i(.hx a x + h z o- l z ) rep- 
resent transitions that do not involve the cavity, and are 
due to the |+) f-> |— ) microwave driving: h x is the mi- 
crowave Rabi frequency, h z is the detuning, and a 1 are 
the Pauli matrices for atom i. In what follows we refer to 
these terms as "fields." Note that the model described 
above, while similar in some ways to the multimode Dicke 
model [35], differs from it in the crucial respect that, in 
the present case, the different modes have distinct spatial 
profiles] it is this feature, not present in the multimode 
Dicke model, that enables frustration to be realized. 

We now proceed to eliminate the cavity modes pertur- 
batively, thus arriving at an effective model for the spins, 
valid on timescales > 2ir/5: 

TT TT , l^( X ')! 2 ga(Xj)g;(x 3 ) i • 

H = F at + 2^ s + ( ' 

a;i<j 

(Note that this result can also be derived from nonequi- 
librium field theory, as in the spinless case [22].) 

Single-mode case. In what follows we denote the ef- 
fective spin-spin coupling as ( = \gfl\ 2 / (A 2 5). Thus, 
e.g., for atoms in a single-mode cavity for which S Q (ir) ~ 
cos(fcx), the zero-field Hamiltonian is 

ffi-modc = (^2. cos(kxi)a\^j(^2 .cos(kXj)aCj+h..c. 

(3) 
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FIG. 2. Phase diagram of frustrated spin systems in cavities, 
as a function of the temperature (vertical axis) and ratio of 
number of modes, p, to number of atoms, N. The thresh- 
olds for self-organization as T — > and as p/N — > were 
computed in Ref. [36], and might have geometry-dependent 
corrections; the boundaries connecting them are schematic. 
Inset: schematic quantum phase diagram for an off-diagonally 
disordered XY model, as a function of hopping (i.e., cavity- 
mediated interaction strength) and spin imbalance, showing 
SF ( "superfluid," i.e., magnetically ordered), BG (Bose glass), 
and RSG (random-singlet glass) phases discussed in the text. 

Because the interaction term can be rewritten as cr^o-j = 
<7 1 x <tI + o- y a y , the system possesses an 0(2) symme- 
try (as the repumping field is negligible). The cavity- 
mediated interaction is ferromagnetic for atoms A apart, 
but antiferromagnetic for atoms A/2 apart; therefore, the 
low-temperature ordered state involves all spins at even 
antinodes aligned along some direction 9 on the equa- 
tor of the Bloch sphere defined by |±), and all atoms 
at the odd antinodes aligned along 9 + it. The interac- 
tions, though disordered (as their magnitude is position- 
dependent), are not frustrated in this case. Note that 
spin-ordering leads to a macroscopic photon population 
in the cavity mode, as in the self-organization of an 
atomic cloud [18, 21, 22, 34] — put differently, magnetism 
is a self-organization of atomic spins rather than posi- 
tions. (This can be seen, e.g., by replacing the a op- 
erators in Eq. (1) by their expectation values.) In the 
driven, dissipative system (with h x ,K ^ 0) such macro- 
scopic occupation corresponds to superradiance. 

Multimode case. We now turn to multimode cavities, 
in which the interactions do not factorize as in i?i- m odc- 
The simplest case is the ring cavity, which supports two 
degenerate modes E±(x) ~ e ±*fc^_ j n this case the inter- 
action term takes the translation-invariant form 

#ring = -C,^COS,[k(x l - X j)] O^O j . (4) 

i<j 

Note that Eq. (4) is precisely Eq. (1) of Ref. [10], which 
approximately describes the RKKY interaction in ma- 
terials (such as Y x Gdi_ x ) having spin susceptibilities 
peaked at a single momentum. While this interaction 
leads to frustration for Ising spins, it does not for XY 



spins; instead, the ground state is a spin spiral of pitch 
A [10, 17]. 

To realize frustration using XY spins, one must 
progress to cavity geometries possessing many degener- 
ate modes, such as confocal and concentric cavities [23]. 
The general Hamiltonian for these is: 

H mm = -( ^ S a(Xi)S a (Xj)cri-(7j. (5) 

For XY spins, Eq. (5) closely resembles the 0{2) general- 
izations [36, 37] of the Hopfield neural-network model [16, 
38]. The mapping to Cook's model [36] is exact for 
translation-invariant, traveling- wave cavity geometries 
such as the ring and confocal cavities; however, the 
basic features of these models (which are similar to 
one another) are expected to extend quite generally to 
Hm-m [10]. 

Associative memories, spin glasses, and self- 
organization. The Hopfield and Cook models describe 
associative memories, consisting of N neurons (i.e., spins 
in the physical system) that collectively encode p "pat- 
terns." In general, p corresponds to the number of cavity 
modes that are resolvable, given the interaction range 
(i.e., x m the notation of Refs. [22]). The associative 
memory is said to function if, starting with any config- 
uration similar to a stored pattern, the dynamics drives 
the configuration to the stored pattern, i.e., if a partially 
self-organized initial configuration at T = becomes 
fully self-organized under the dynamics (this point is 
discussed further in Ref. [39]). In the Hopfield and 
Cook models [16, 36, 38], this is the case (as TV — > oo) 
for small p/N(< 0.05), e.g., in the single-mode cavity. 
For p/N > 0.05, metastable states proliferate, and the 
system becomes a spin glass; the spin glass differs from 
the self-organized phase in that the ground-state atomic 
configuration does not globally emit superradiantly 
into any particular cavity mode; nevertheless, it is a 
distinct phase from the high-temperature paramagnetic 
phase [9] . A finite-temperature phase transition between 
the two is known to exist in the case that p/N — > oo, 
i.e., the Shcrrington-Kirkpatrick model [9, 16]. These 
considerations lead to the global phase diagram shown 
in Fig. 2. 

Tuning and detection. Both the associative memory 
and the spin glass are low-temperature phases. The for- 
mer is stable when fcgT < h(N. As is standard in ultra- 
cold atom experiments, the temperature is determined 
(in the n <C S, 7 <C A limit) by the system's initial en- 
tropy; however, £ increases with pump laser intensity, 
and can be tuned across the transition. The spin order- 
ing threshold is similar to that for self-organization, and 
is achievable, even for relatively large S, for reasonably 
long experimental lifetimes [20-22]. The effective num- 
ber of modes coupling to the atoms can be decreased by 
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adjusting the length of the cavity away from the confo- 
cal/concentric limit; as £ ~ 1/5, only modes having suf- 
ficiently small 5 couple strongly to the atoms. The spin 
glass transition temperature T g in Hopfield-type mod- 
els is comparable to the single-mode ordering tempera- 
ture [16]. 

The self-organized phase should be detectable via the 
light emitted from the cavity, but the spin-glass phase is 
not, as it does not exhibit superradiance. One straight- 
forward way to detect this phase is through its slow relax- 
ational dynamics: a possible protocol involves initializing 
all spins in a certain region in the |+) state via a lo- 
cal spin addressing protocol [40] , and measuring the spin 
relaxation timescale (observed, e.g., via phase-contrast 
imaging [41]) as a function of pump intensity. A fea- 
ture common to both the superradiant and spin-glass 
phases at low temperatures is the presence of a large num- 
ber of low-energy excitations; these reveal themselves in 
condensed-matter systems via the heat capacity. In the 
cavity QED setting, such excitations can be detected, 
e.g., via two-photon spectroscopy [42]. Further possibil- 
ities for distinguishing the two low-temperature phases 
via their response functions are considered in Ref. [43]. 

Quantum regime. Thus far, we have focused on the 
classical spin physics realizable using cavity-mediated in- 
teractions. We now turn to the quantum regime, in which 
Eq. (5) can be mapped [44] onto a Bose-Hubbard model 
in the limit U/t — > oo, via the transformation er + — > &L 

Hbh = -wY^t ij (b\b j +h.c.)+iM^2b\b i . (6) 

ij i 

According to this mapping, a | +) state corresponds to the 
presence of a b boson whereas a |— ) state corresponds to 
the absence of a b boson; the chemical potential n is de- 
termined in the standard way from the number of bosons 
(i.e., |+) atoms). More generally an n-state atom maps 
onto a Bose-Hubbard model with a maximum occupation 
per site of n — 1: 



interact identically with the trapping lasers). This is 
challenging to achieve in optical lattices [26], as varying 
the hopping amplitude via the lattice depth inevitably 
leads to on-site disorder. The phase structure of the dis- 
ordered Bose-Hubbard model is known to be richer in 
the absence of chemical potential disorder, especially in 
one dimension [24]: the off-diagonally disordered model 
exhibits a Mott glass phase, as well as a random-singlet 
glass phase [27, 28] (see Fig. 2). These phases arc not 
present in models having on-site disorder, and are thus 
not directly realizable in optical lattices, but are realiz- 
able in the cavity-based setting. 

A simple one-dimensional geometry that realizes the 
model of Ref. [24] involves a chain of atoms trapped 
perpendicular to the cavity axis. The cavity modes are 
Hcrmite-Gaussian along this direction, which we label 
y; thus, the nth mode has a profile of the Hermite- 
Gaussian H n (y/L) exp(—y 2 /L 2 ), where L is the waist of 
the TEMoo mode. For good confocal cavities n ~ 10 — 
100; a more scalable geometry involves atoms trapped 
along the cavity axis. An atom at position y couples 
most strongly to modes with "classical turning points" 
near y; these modes have large amplitudes within a dis- 
tance L from the turning point, and either decay or oscil- 
late rapidly beyond this distance. Thus the interaction 
range is given by L and (in particular) is finite, and the 
analysis of Ref. [24] applies. For the case of a A atom hav- 
ing equal spin populations, one can realize the random- 
singlet glass, which possesses long-range spin-singlet cor- 
relations [27, 28]. (In order to realize the Mott glass, one 
could use, e.g., an atom possessing three ground states.) 
The three realizable glassy phases can be distinguished, 
e.g., via internal-state dependent transport or compress- 
ibility measurements [26]. 
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ffn-icvci = -w^tij{b\bj + h.c.) + U ^(n t - n) 2 , (7) 

ij i 

where the "interaction" term U can arise, e.g., because 
of the quadratic Zeeman shift. Note that, in this bosonic 
terminology, a "superfluid" state corresponds to a finite 
expectation value of o± (i.e., to in-plane magnetic order- 
ing). Thus, it is a "superfluid" of spins and not of the 
atoms, which are frozen in place. 

A crucial difference between cavity-based realizations 
of the Bose-Hubbard model and optical-lattice ones [26] 
is that, in the cavity-based setting, strongly disordered 
hopping amplitudes are natural (hopping amplitudes be- 
ing determined by the; oscillatory cavity mock; functions) 
even in the absence of on-site disorder (as both spin states 
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